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Abstract. We consider the spherical model on a spider-web graph. This graph is 
effectively infinite-dimensional, similar to the Bethe lattice, but has loops. We show 
that these lead to non-trivial corrections to the simple mean-field behavior. We first 
" determine all normal modes of the coupled springs problem on this graph, using its large 

symmetry group. In the thermodynamic limit, the spectrum is a set of ^-functions, 
and all the modes are localized. The fractional number of modes with frequency less 
(— i ■ than u varies as exp(— C /uj) for uj tending to zero, where C is a constant. For an 

O . unbiased random walk on the vertices of this graph, this implies that the probability 

^ , ' of return to the origin at time t varies as exp(— C't 1 / 3 ), for large t, where C is a 

constant. For the spherical model, we show that while the critical exponents take the 
values expected from the mean-field theory, the free-energy per site at temperature T, 
near and above the critical temperature T c , also has an essential singularity of the type 
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1 Introduction 

It is well-known that in systems having dimension d greater than the upper critical 
dimension djj, fluctuations in thermal equilibrium near critical phase transitions are 
well described by the mean-field theory pQ . For simple Ising-like spin models, the value 
of upper critical dimension is known to be 4. In this paper, we give an example of 
a lattice, whose effective dimension is infinite, and the exact solution of the spherical 
model on the lattice shows that while the critical exponents are correctly given by 
the mean-field theory, (e.g. the specific heat has a discontinuity), the free-energy near 
the transition also shows an essential singularity, with the singular part varying as 
exp[— K{T — T c ) -1 ^ 2 ], for temperatures T near and greater than the critical temperature 
T c . While it is known that higher derivatives of the free energy do show logarithmic 
singularities in the spherical model for even dimensions d > du [2], this seems to be the 
first case where an essential singularity is encountered in the sub-leading singular part 
of the free energy in an infinite dimensional model. 

The lattice we study is the spider-web lattice, also known as spider-web graph or 
the spider-web network in literature. These graphs were introduced by Ikeno in 1959, 
for efficient design of crossbar exchanges for telephone switching networks [3]. The name 
spider-webs to describe them seems to have been first used by Feiner and Kappel in 1970 
[3]. With the advent of digital switching and mobile communication networks, crossbar 
exchanges have become obsolete, and their optimal design is no longer of technological 
interest. But the spider-web graphs are very interesting structures for statistical physics. 

A 2 N x M - node spider network graph consists of M-levels each of which has 2 N 
vertices. Each vertex in the m-th level is connected to two vertices each in the (m± l)th 
levels. These are finite graphs where all vertices are equivalent. In the thermodynamic 
limit when the number of vertices in the graph tends to infinity, the number of sites 
within a distance R of a given site increases exponentially with R, and the effective 
dimension of these graphs is infinite. These are thus similar to the Bethe lattice, which 
has been studied a lot in statistical physics. However, for the Bethe lattice, realized 
as the limit of tree graphs, separating the bulk properties from the surface properties 
requires care as most of its sites are near the surface. One can realize a Bethe lattice 
without surface effects as the limit of a random graph with uniform coordination number, 
but this involves averaging over disorder, which may be quenched [5], or annealed [6]. 
Spider-web graphs avoid the difficulties associated with such averagings. They also have 
additional symmetries, in addition to translational invariance, which makes the study 
of statistical physics models on these graphs analytically more tractable. 

Given its origins, it is not surprising that most of the earlier studies dealt with 
the percolation properties of these networks [3]. Takagi |7] showed that in a certain 
class of graphs, the spider-web network has the highest linking probability. But, Chung 
and Hwang [8] showed that, outside this class, there are other graphs that are better. 
For subsequent work on the percolation properties on these networks, see [HI HOI EEU H2j • 
Other models on these graphs have not been studied much so far. There is large amount 
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work dealing with spectra of different graphs. Some basic references may be found in 

In this paper, we study the spherical model on the spider-web graph. We start 
by determining the normal modes of vibration of masses coupled by nearest -neighbor 
harmonic springs on this network. We obtain the entire spectrum for finite values of M 
and N, and also describe the qualitative features of the spectrum in the thermodynamic 
limit of M — > oo, with N finite. In the double limit of both M, N — > oo, almost all the 
eigenmodes of the springs problem are localized on this network. We study the return 
to origin for the related problem of a random walker on this network. We show that the 
probability that a random walker returns to origin after t steps varies as exp(— 1 1//3 ). This 
should be compared with a power-law decay t~ d / 2 on d- dimensional hypercubic lattices, 
and an exponential decay on the Bethe lattice. We use the knowledge of normal modes 
of the springs model on this lattice to determine the critical behavior of the spherical 
model on this graph, and show that the free energy shows an essential singularity at the 
transition point. 

The plan of this paper is as follows. In section II, we define spider-web graphs, and 
then define the springs Hamiltonian, the random walk problem and the spherical model 
on these graphs. In section III, we discuss the decomposition of the Hilbert space of 
the displacement field in the springs problem into subspaces, called sectors, using the 
symmetries of the spider-web graphs. In Section IV, we determine the frequencies of all 
the normal modes of the springs Hamiltonian, and also determine the eigenmodes. In 
Section V, we use this spectrum to determine the exact generating function for returns 
to origin of a walker doing a random walk on the vertices of the spider web graph. In 
section VI, we show that the singular nature of density of states of the springs model for 
small frequencies implies that free energy of the spherical model shows a discontinuity 
in specific heat, as expected from the standard mean-field theory, but also an essential 
singularity. 

2 Definition of the model 

A spider-web network is a graph consisting of Mq N vertices, divided into M levels, each 
level with q N vertices, with M, N and q being positive integers. We label the levels by 
integers 1 to M, and vertices within a level are labeled by iV-digit integers r, base q, 
with < r < (q N — 1). In the following we discuss only the case q = 2. Generalization 
to other values of q is straight-forward. 

On a binary string of length N, we define an operation R of adding a zero from 
the right, making space for it by shifting all bits by one space to the left. The left-most 
(first) bit of the original string is lost. So, for example, for N = 6, action of Rq on 
the binary string 111001 yields a string 110010, and we write i? (H1001) = 110010. 
The operator Ri adds a 1 to the right. We similarly define operators Lq and L\, which 
correspond to adding and 1 respectively at the left end, and shifting the bits one space 
to the right. In this case, the rightmost bit is lost. Clearly, Ro{r) = 2r(mod 2 N ), and 
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#x(r) = 2r + l(mod 2 N ). Also, L (r) = [r /2J and Li(r) = 2 JV_1 + [r/2j. 

Now consider a vertex labeled by a binary integer r in the level m. We denote this 
vertex by (m, r). The spider-web graph is defined by connecting this vertex to the two 
vertices (m + 1,Rq(t)) and (m + l,Ri(r)) in the level (m + 1). We assume periodic 
boundary conditions, and the vertices in the level m = M are connected to level m = 1. 
It is convenient to assume that the index m is defined only modulo M, so that m — 1 is 
the same as m — M + 1. We will take these connecting edges to be undirected. Then, 
each vertex (m,r) is also connected to the vertices (m — l,Lo(r)) and (m — l,Li(r)). 
We will assume that M > N + 1, as otherwise the graph is not fully connected. Fig. [1] 
shows a graph with M = 4, iV = 3. 

Given this structure of connections, it is easy to see that there is unique directed 
path of length TV from any vertex (m, r) to any of the 2 N vertices in the level (m + N). If 
the binary label of the target vertex is iii 2 i3-..iN, one just starts at the vertex (m, r), and 
takes steps to the next level using edges corresponding to R il , R i2 , R i3 . . . respectively. 

In the context of crossbar exchanges, we think of each node of the graph as 2 x 2 
crossbar, with two input lines, and two output lines. It is natural to work with open 
boundary conditions, with no direct connections between level M and level 1. Then 
there are 2 N+l input lines coming into level 1, and there is an equal number of lines 
coming out of level M. The total numer of crosspoints used is 4 per vertex of graph. 
Then, using only M = N levels, we can connect any of 2 N+1 input lines into the level 
m — 1, we can connect to any of the 2 N+1 lines of the crossbars in the level m = N. 
Also, there is a simple algorithm to determine the path from the input node to the 
output node. This is much less than the number 2 N+1 x 2 N+l , which would have been 
required, if a single stage crossbar was used. If M > N, we get multiple possible paths 
from the input to output. This was the original reason for interest in these networks. 

For studies in statistical physics, one is usually interested in the thermodynamical 
limit, where the effect of boundaries can be ignored. We note that here, local properties 
of the network do not depend on M or N, for distances up to N. For example, the 
number of sites that are within a distance r of any site is the same for all sites, and is 
independent of M and N, so long as N > r. Thus, for systems with interactions of finite 
range, say for a nearest-neighbor Ising model defined on this graph, the thermodynamical 
limit will be well-defined. As another example, if there is a random walker that moves 
on the vertices of this graph, taking a step at random along one of the four edges at the 
current position, the probability of return to starting point after T steps is well defined, 
and independent of M and iV , so long as T < 2N. 

We consider a unit mass attached to each vertex of the spider-web network, coupled 
to its four neighbors by harmonic springs of unit spring constant. We denote a scalar 
variable x(v) for each vertex v = (m, r) of the graph. Let the corresponding momentum 
variable be denoted by p(v). Then, the classical -mechanical springs Hamiltonian is 
given by 

H = \ 5>( v ) 2 + \ E *( V ) K ( V > v >( v ') (!) 

v v,v' 
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Figure 1. N = 3, M = 4 spider network with periodic boundary conditions, the 
vertices in the column m — 5 arc identified with those of m = 1. A vertex with 3-bit 
label £ in the level m is connected to the vertices with labels Rq(£), and R\{1) in the 
level (m + 1). 



where K is the matrix of coupling constants defined by 
K(v, v') = 4, if v = v', 

= — 1, if v, v' are nearest neighbors, 

= 0, otherwise (2) 

Clearly, the matrix IK can be expressed in terms of the adjacency matrix of the 
graph. As is well-known, the problem of determining the eigenmodes of the set of coupled 
oscillators reduces to diagonalizing the interaction matrix K, equivalently finding the 
spectrum of the discrete Laplacian of the graph. Let Va(v) be a normal mode of this 
system with frequency oj a . Then, ip a is a solution of the eigenvalue equation 

K^ a (v)=w^a(v). (3) 

We will denote by F(u 2 ), the fractional number of modes of an eigenvalue less than 
or equal to u 2 . Hence, < F{u 2 ) < 1. This is called the cumulative spectral density. 
The spectral density D(lj 2 ) is defined as dF(cu 2 )/duj 2 . In the limit M —> oo, with N 
finite, we will show that D(co 2 ) has both a discrete part, and a continuous part. The 
fraction of modes that give the continuous part is 2~ N . Therefore, in the limit N — > oo, 
we get only a discrete spectrum. 

It is well-known that the diagonalization of the springs Hamiltonian is also related 
to problem of random walks. Consider a continuous-time random walk on the vertices of 
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this graph, where the walker can jump to any one of the neighboring vertices with rate 
1. Thus, in a short time dt, it can jump to any of the four neighbors with probability 
dt, and stay at the same vertex with probability 1 — Adt. Suppose the walker starts at 
vertex v . Then, we denote the probability that the walker is at vertex v at time t, by 
Prob(v, t). We define a vector \P(t)), whose v-th component is Prob(v, t). Then the 
evolution of \P(t)) satisfies the Markov equation 

j t \P(t)} = -K\P(t)} (4) 
Then, the probability of being at vertex v' at time t is given by 

Prob(v',t|vo,0) = (v'|e" K4 |vo). (5) 

The probability of being at the starting point at time t is given by Prob(vo, t|vo, 0). 
From the translational symmetry of the graph, this does not depend on vo- We define 
the Laplace -transform of return probabilities R(s) by 

POO 

R(s)= dtProb(v ,t\v ,0)e- st (6) 
Jo 

Then it easy to see that 

1 r°° 

R{s) = m zTr[a + K\- 1 = J o du 2 D(u 2 )[s + c 2 ]- 1 (7) 

We use this determination of the normal modes to solve the spherical model on 
this graph. In the spherical model, we assign a scalar variable s v to each vertex v of 
the graph. These variables are subjected to a single constraint, known as the spherical 
constraint 

£4 = #. (8) 

V 

where N s is the total number of sites in the graph. The spherical model can be obtained 
as the large n limit of an n- vector model on the same graph [2], where we assign an 
n- vector spin S v to each site v, with a constraint |5 V | 2 = 1, at each vertex v. The 
hamiltonian of the spherical model is given by 

v,v' 



3 Symmetries and sector decomposition 

As noted before, the problem of determination of the normal modes of the springs 
hamiltonian H reduces to that of diagonalizing the matrix IK. This is an N s x N s 
matrix, where N s is the number of vertices in the graph, N s = M2 N . The matrix K 
may be viewed as the discrete laplacian operator, acting on the vector space of square- 
integrable real functions defined on the vertices of the graph. Using the symmetries 
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of the problem, we can block-diagonalize the matrix K. We will call the subspaces 
corresponding to different blocks as sectors. 

The spider-web network has a high degree of symmetry. Firstly, all vertices are 
equivalent. We define a translation operator T which takes a vertex (m, r) to (m + 1, r). 
Clearly, T is a symmetry of the graph, and we have 

T M = I (10) 

where I is the identity operator. 

There are additional symmetries. Consider any integer % lying between 1 and 
M. Define fk{ r ) as the integer obtained from r by flipping the fc-th bit of its binary 
representation, reading from left to right. If the bit is 0, it is replaced by 1, and vice 
versa. We define the operators F&, for k — 1 to N. The operator F& takes the site 
(m + k — l,r) to another site (m + k — 1, fk{f)), for all r,m. Sites in levels j where 
m — j + 1 (mod M) is outside the range 1 to N are left unchanged. Clearly, we have 
F^ = I. Also, we see that this mapping preserves the adjacency strucure of the network: 
if vertices v and v' are nearest neighbors, so are F,(v) and Fj(v'). Hence we have 

pK,Fi] = 0. (11) 
Acting on the space of wavef unctions, the action of the operators Fj is given by 

F^(v) = ^(F,v) (12) 
Now, clearly, the operators T and Fj, % = 1 to M commute with "K: 

pK,Fi] = [K,T] =0. (13) 
Also, different Fj commute with each other: 

[FijFj-] = 0. for all i,j. (14) 

We define shift-operator M. , which acts on the set of vertices of the graph, and for 
all sites (m,r), M.o(m,r) = (m + l,i?o( r ))- Operators Ki,Lo,Li are defined similarly. 
Then clearly, we have 

K = L + Li + R + Mi - 4 (15) 

We can think of Lo, Li, Mo and Mi as the generators a semigroup. These can be written 
as iV s x N s matrices. 

If we first apply E a , and then on any vertex v, where a and b take values or 1, 
the result is independent of a. This implies that Mo^o = M.oI j i- More generally, we have 

IRo-L'O = MqILi, K.1L0 — IR1IL1, ILqMq = ■LoMi, IL^Mq = IL1IR1 (16) 

Similarly, it is easily seen that for all allowed values of the subcscipts a, b, c 

R c L, b R a = K c , lL c R b h a = L c . (17) 

One can write similar equations involving longer products of the type IRIRLLIRIR etc.. 
This is interesting for an purely algebraic characterization of the spider-web graphs, but 
it is not needed for our discussion here. 
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We can find eigenvectors tjj that are simultaneous eigenvectors of K and Fj, for all 
i, 1 < i < M. 

Now, we denote the eigenvalues of Fj by /j. These are ±1. Consider an eigenfunction 
i/j, with 

Wiip = frf, for i = 1 to M. (18) 

The Hilbert space of all functions ip(m,r) breaks into 2 M disjoint sectors, 
each corresponding to a subspace spanned by basis vectors that are eigenvectors 
corresponding to a specified set {fi}. For M 3> N, the number of sectors can be 
substantially greater than the M2 N , which is the total dimension of the Hilbert space. 
It follows that for M ^> N, a large number of sectors must be null. 

It is easy to see that if ip(m,r) is an eigenfunction of Fj with eigenvalue — 1, 
then if;(m,r) must be zero, whenever m [i — N + l,i}. If more than one of the 
/j's are negative, the corresponding eigenfunction ip(m,r) can be nonzero only on the 
intersection of the corresponding intervals. This severely constrains the allowed values of 
{fi} for a non-null sector. In particular, this implies that any sector with Fj = Fj = — 1, 
and distance between i and j greater than N — 1 must be null. 

Only in the sector with all /, = +1, the sector decomposition does not break 
translational invariance, (the sector with all f = — 1 is zero- dimensional for M > N). 
In this case, we can break the sector further into M one-dimensional subsectors using 
the eigenvalues of T. 

It is straight forward to work out the sector decomposition and the spectrum, for 
small values of N explicitly. 

3.1 The case N = l 

For N = 1, the graph is shown in Fig. [2j There is one sector with all fi = +1. This must 
have ip(m, 0) = ip(m, 1). Thus there are M linearly independent variables in this sector, 
and dimensionality of this sector is M. These are further divided into one- dimensional 
subsectors using the eigenvalues of T. 

Then, there are M sectors with exactly one f = —1, and rest positive. In the 
sector with fi = —1, we must have ^(1,0) = —^(1,1), and i[>(j,0) = = for 

j^l. So, each of these sectors is one- dimensional. The corresponding eigenfunction is 
fully localized within one level, and is shown in Fig. |2j All sectors with more than one 
fi negative are empty. The total dimension of the space is thus = 2M, equal to total 
number of modes in the problem. 

3.2 The case N = 2 

For N = 2, the sector with all /j's positive has dimension M. It is divided into M 
one-dimensional subsectors using the eigenvalues of T. 

There are M sectors where only one / is negative. Consider the sector with f\ = — 1, 
and all other /'s positive. In this case, let the eigenfunction be ip(m,r). We have 




Figure 2. Figure shows a part of the N = 1 network, with a localized eigenfunction 
corresponding to exactly one of the Fj's negative. 



¥iip(m,r) = —ip(m,r). If m is not equal to or 1, we have ¥i(m,r) = (m, r). 
This implies that ip(m,r) must be zero, for m ^ 0,1. Hence the only non-zero 
elements are ^(0, 0), 0(0, 1), ^(0, 2), 0(0, 3) and V(l, 0), 0>(1, 1), ^(1, 2), 70(1, 3) . Using 
the symmetries of Fm and F 2 we also have -0(0, 0) = -0(0, 2) = —"0(0, 1) = —"0(0, 3) and 

0) = 1) = —0(1, 2) = — 0(1, 3). Thus there are only two independent elements 
and hence this sector is of dimension 2. 

The number of sectors with exactly one f^ a negative is clearly M. Each of these is 
two dimensional. 

Now consider the sectors with exactly two of the /'s negative. Say, we have 
fi = fj = —1, and the rest of the / are positive. Clearly, we can assume % ^ j. 
Using the result proved above, if \z — j\ > 1, then the sector is empty. 

Hence, the non-trivial cases with two of the f's negative are only with fi = f i+ i = 
— 1, for some %. Clearly, there are M choices of i. Say % — 1. then the corresponding 
eigenfunctions can be non-zero only for sites in the level m — 1. The values of 
ip{m = l,r) for different r can be related to each other using the the symmetries of 
Fi and F 2 . Thus this sector is one-dimensional. 

If at least three of the /'s are negative, then one can find a pair of values 
fi = fj = — 1 , with \i — j\ > 1, and the sector must be empty. 

Thus, we find that for N = 2, the only non-empty sectors are those with 
total number of negative eigenvalues of Fj at most 2. The total number of sectors 
M + 2M + M = AM, which is equal to the number of different vertices in the graph. 

3.3 Higher N 

It is straightforward to extend the treatment given above to higher N. We will assume 
M 3> N, for convenience. 

Given all the eigenvalues {fi}, we can relate the value of the eigenfunction ip(m, r) 
for any r to the value at the vertex (m, 0) by applying the operators Fj. Clearly, it differs 
from ip(m, 0) by at most a sign, and ijj(m, r) = ±<f) m . The dimension of the sector is the 
number of nonzero m 's that are allowed. 

Consider first the sector with all fi — +1. The corresponding eigenfunction ip(m, r) 
must be of the form ip(m,r) = g m , independent of r. There are M possible choices of 
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m. Hence this sector is M- dimensional. In this sector, T is a good symmetry, and we 
can use the the eigenvalues of T = tj = exp(i-?f), with j — 1 to M, to further divide 
this sector into M one-dimensional sectors. 

Now consider the sectors with exactly one fi negative. Say fi = —1. Then the only 
nonzero values of 4> m are for m = 1, 0, —1, —2, . . . , 2 — N. This sector is of dimension 
N. There are M such sectors. 

For more than one negative fi's, if the largest distance between two negative /'s 
is k, then k < N — 1. The number of non-zero g m 's can be only N — k. Hence this 
is the dimension of each such sector. There are k — 1 /'s between the outermost two 
negative /'s in this band, which can be ±1, and so the number of sectors with a given 
position of leftmost and rightmost negative /'s is 2 k ~ l . also the left most negative / has 
M possible positions. Hence the number of sectors with a given value of k is M2 k ~ 1 , 
and the dimension of each sector is N — k. 

As a check, this gives the dimension of the full Hilbert space of eigenfunctions as 

AT-l 

M + MN + M^2(N- A;)2 fc - 1 = M2 N (19) 

k=l 

4 Spectrum of normal modes of the springs hamiltonian 

The problem of diagonalization of K is considerably simplified using its block- 
diagonalization induced by the symmetries {Fj,T}. It is instructive to consider the 
cases of small iV first. 

4-1 The spectrum for N =1 

In this case, K is fully diagonalized using the symmetries alone. There are 2M 
eigenstates, and the 2M nontrivial sectors are all one-dimensional. We can thus write 
down the eigenvectors by use of symmetries alone, and then determine their eigenvalues. 

There are M states in the sector with all j\ = +1. These are further labeled by 
eigenvalues of T. The vector corresponding to the eigenvalue e lk of T is given by 

ip(m, 0) = ip(m, 1) = e ikm , (20) 

with k = (27rj)/M, and j taking integer values to M — 1. The corresponding 
eigenvalue of K is easily seen to be 

uj 2 = A-Acosk (21) 

Other non-null sectors have exactly one fi = —1, and rest positive. The eigenvector 
corresponding to the f] = — 1 is easily written down 



^(i,0) = -^(i,l) = l; (22) 
i/j(m, r) = 0, for m ^ i. 
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This mode is fully localized on the sites in the level i. The corresponding eigenvalue 
of K is seen to be 4. 

In the limit M — > oo, we see that the localized modes are degenerate, and give a 5- 
function peak in the spectrum D{uj 2 ) at u 2 = 4. The total weight of this peak is 1/2. 
The remaining modes are seen to be extended, and give rise to a continuous spectrum. 
Adding the delta-function peak to the continuous part given by (l2*Tj) . the total density 
of states is given by 

D (w ») = fa - 4) + 4tK(8 1 _^ )|1/2 ; for N = 1 (23) 
4-2 The spectrum for N = 2 

For the sector with /j = +1 for all i, again we have M modes. These are given by 

ip(m,r) = exp(ikm), independent of r, (24) 

and the corresponding k values , and the eigenvalue of K is as given by (121j) . 

Now consider the sector with only one negative /, say fa = —1, and the rest of 
/'s are all positive. This sector is two-dimensional. From the symmetries alone, we see 
that the corresponding eigenfunction ip(m, r) is non-zero only for m = i and m — i — 1. 
Using /■_! = f i+1 = +1, we get 

^(i, 00) = 01) = -W, 10) = -^(z, 11) = ^(i) (25) 

- 1, 00) = -V>(« - 1, 01) = ^(^ - 1, 10) = -V>(i -1,11)= ^(i - 1) 

The eigenvalue equation relating p(« — 1) and is easily seen to simplify to 

-u 2 g(i-l) = -4g(i-l) + 2g(i), (26) 

-u 2 g(i) = -Ag{i) + 2^(i - 1) 

This is a coupled set of two linear equations. It is easily seen that allowed solutions 
are g(i) = ±g(i — 1), and the corresponding values of u 2 are 2 and 6. The eigenmode 
corresponding to u 2 = 2 is shown in Fig. [3j 

For the sector with fi = = —1, there is only one eigenvector, which is fully 
localized in the level m = i. The corresponding eigenvalue is 4, and the eigenvector is 
given by 

^(;, 00) = -i>(i, 01) = 10) = ^(i, 11) (27) 
In the limit M — > oo, the density of states is seen to be 

D{uj 2 ) = \[5{uj 2 -2)+5{u 2 -A)+5{u 2 -Q)] + ^ \ ■ for N = 2(28) 




Figure 4. An cigenmode in the network with N > 2 with eigenvalue 4. The 
amplitude, shown in red, is nonzero only at four sites 



4-3 The spectrum for higher N 

The treatment given above is easily extended to higher N. Firstly, given the set of 
values {fi}, the corresponding eigenfunction ip(m,r) is of the form ip(m,r) = ±g(m), 
and the sign is easily determined in terms of {fi} and the binary representation of r. 
We choose the convention that the sign is positive for r = 0. 

A nontrivial eigenfunction is possible, only if the indices % corresponding to negative 
/'s occur in a narrow band of width at most N. If the width of the band is £, the number 
of non-zero g(m) is u = N — t + 1. 

Consider, for simplicity, the sector in which the non-zero g(m) y s are for m G 
[i, i + u — 1]. It is easily seen that the eigenvalue equation (EJ) in terms of the function 
g(m) becomes 



- cu 2 g(j) = 2g{j + 1) + 2g(j - 1) - Ag(j), for j e [i, z + u - 1], (29) 
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= 0, otherwise. 
This equation is easily solved. The eigenvalues are 

ITS 

w 2 = 4-4cos( ), for s = 1 to m. (30) 

u + 1 

and the corresponding eigenfunction is 

/ n a ■ ^(m - i + l)s . . . 

g 8 {m) = Asm[ — ^ J, lor m G [z, i + it — 1J, (31) 

= 0, otherwise. 

As a simple example, if fi = J'n = — 1, and all /j = +1, for j > N, then the 
eigenfunction g(m) is non-zero only for m = 1, and the eigenvalue is w 2 = 4. There are 
2 7V_2 such sectors, corresponding to different possible choices of fi for i G [2, iV — 1]. All 
these different modes are localized in the level m — 1. It is then possible to construct 
linear combinations of these modes, which are still eigenfunctions of K, but no longer 
eigenf unctions of Fj. We find that these can be chosen so that each eigenfunction is 
highly localized, and nonzero only at four sites. 

This mode is shown in Fig. HI Let S be a binary string of length N — 2. We consider 
the eigenfunction which is nonzero only at the four sites (1, 050), (1, 150), (1, 051) and 
(1, 151), defined by 

ip(l, 050) = -^(1, 051) = -^(1, 150) = ip(l, 151) (32) 

Clearly, this is an eigenfunction of K with eigenvalue 4. 

A similar construction works for all other localized eigenmodes. An eigenmode 
which extends over u levels will have a degeneracy of 2 N ~ U ~ 1 , and one can choose a 
linear combination that localizes the mode maximally. Such a mode will have non-zero 
amplitude only at 2 n+1 sites in each of the u layers. 

It is easily seen that the spectral density function for a general iV is: 

N-l 1 u 

f(- 2 ) = E^£^ 2 -( 4 - 4 <-(^)>) 

u=l s=l 

1 N 

(33) 

In the limit of M, N — > oo, the spectrum becomes purely discrete. This is very 
interesting, as the only other known example of a regular transitive infinite graph with 
a discrete spectrum of the laplacian is the Cayley graph of the lamplighter group, or its 
generalizations 
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Figure 5. The spectral density in the limit of large N. Only the peaks with weight 
greater than 10 have been shown. 



5 Returns to the origin 

From the expression for D(u 2 ), it is straightforward to determine how the probability 
of the random walker being found at the starting site at time t varies with t for large 
t. Using (jHJ), and the fact that the return probability Prob(v , t; v , 0) does not depend 
on vo, it may be written as 

POO 

Prob(v ,t; v ,0) = / du 2 D(uj 2 ) exp(-tuH) (34) 
Jo 

For finite N, with M large, the network is effectively one- dimensional. The 
integral above is dominated by by the continuous spectral density for small u 2 . As 
D{oj 2 ) ~ 2~ N ^ for small u, it is seen that Prob(v ,t; v , 0) varies as 2~ N t~ 1 / 2 for large 
t. 

For N large, the spectrum is purely discrete. In ( 130]) . For a given value of the 
localization width u, the mode with lowest frequency corresponds to r = 1, and then 
for large u, the frequency may be approximated by u 2 ~ 2ir 2 u~ 2 . Then, we can write 

Prob(v ,t;v ,0) ~ V2-"exp(-^) (35) 

u 

For large t, the leading behavior of this summation can be determined by steepest 

1/3 

, and we get 



descent. The term corresponding to u = u 

3 



In 2 



Prob(v , t; v , 0) ~ exp 



(4^)1/3^1/3 



for large t. (36) 



We note that this is qualitatively different from the behavior for the Bethe lattice, 
where this probability decays exponentially with time. 

The function R(s) in the limit iV — > oo, using (JZJ) and ( 133|) . is given by the formula 

oo u ^ 

^)-EE 2 , + i (g + 4-4cos^ T ) (37) 

u=l r=l K w+1' 
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We note that R(s) can be written as 

R{s) = (l/2)Q(s) - (38) 

where 

«»> = E E 2^.(, + 4 -4c M -^) (39) 

u=l r=l y 

The summation over r in the above equation can be done explicitly using the 
identity 



2u+l 

y — - 

^— ' 1 — a; cc 



2(n + 1) 



« +1 + a u _ +1 ) 
{a u + +1 - a u _ +1 ) 



(40) 



T n+1 (l/x) 

U n (l/x) { ' 

where a± = (l±vl — x 2 )/x and T n and U n are Tchebychef polynomials of the first 
and second kind respectively [T5j . 

6 Critical behavior of the spherical model 

The calculation of equilibrium properties of the spherical model is quite standard [2]. 
We take care of the spherical constraint by using a Lagrage multiplier, and consider the 
Hamiltonian 

n sph (\) = H sph + ^\J2 s l ( 42 ) 

V 

where A is temperature-dependent constant, whose value is determined by the condition 
5>v>A = iV,. (43) 

V 

Here ( )> denotes equilibrium thermodynamical average using the Boltzmann factor 
corresponding to the Hamiltonian H sph (\). Since T-L sp h(X) is a purely quadratic 
hamiltonian, in terms of normal mode oordinates, it becomes uncoupled oscillators. 
Thus, we get 

a " 

It is easy to see that A is a monotonic non-decreaseing function of T. As the 
temperature is decreased from very high values, at a particular temperature T c , in the 
thermodynamic limit A becomes zero, and sticks to the value zero for all T < T c . The 
transition temperature T c is given by the equation 



oo 

k B T c I du 2 D(uj 2 )/cu 2 = 1 (45) 



o 
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For T < T c , the spherical constraint Eq. (jSJ), is satisfied by by a macroscopic 
occupation of the uniform mode. If the uniform magnetization is m(T), it is easy to see 
that we have 

m 2 (T) = 1 - T/T c , for all T < T c . (46) 

For T near T c , this implies that m(T) ~ (T c — T) 1 / 2 , and the magnetization exponent 
f3 takes the classical value 1/2. The internal energy per mode is the classical value |/csT, 
by the equipartition theorem, and the specific heat takes a constant value (1/2) ks for 
all T < T c . 

For T/T c = 1 + e, with e small and positive, it is easy to see that A is a linear 
function of e: 

A = Ce + £>(e 2 ) (47) 

where C is a positive constant. The specific heat for e > is smaller than /c_b/2, and 
jumps discontinuously at T c , with the discontinuity proportional to the C. 
The free energy per site f(T) for T > T c is given by the equation 

f(T) = f du*D^) MttTt) (48) 

2W S J ur + A 

Using the fact that D(u 2 ) varies as exp(—C'/u) for u tending to zero, we see that 
for small positive e, f(T) can be written as sum of a Taylor expandable regular part, 
and a singular part f S in g (T) ~ exp(— C'/y/e). 

It would be interesting to extend this analysis to other models of phase transitions 
on this lattice. It seems reasonable to expect somewhat similar behavior. For example, in 
the n-vector models with nearest -neighbor ferromagnet couplings, the high-temperature 
expansion for free energy involves sum over loop diagrams. The number of self-avoiding 
polygons of a given perimeter i may be expected to have a behavior similar to number 
of random walks that return to origin. Note that these essential singularities can not be 
captured easily in a perturbative approach. Another interesting question is the effect of 
disorder on the localization properties of the normal modes of the springs Hamiltonian 
on this network. These seem to be interesting directions for further study. 
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